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ABSTRACT 

Two  aspects  of  Incompressible  laminar  boundary  layer  flow  on  a  semi- 
infinite  flat  plate  at  zero  incidence  are  discussed.  Firstly,  Alden's 
proposed  scheme  of  solution  of  the  Navier-Stokes  equations  for  this 
problem  in  terms  of  an  expansion  in  inverse  powers  of  the  local  Reynolds 
number  is  demonstrated  to  be  invalid,  inasmuch  as  the  boundai'y  conditions 
in  the  free  stream  cannot  be  satisfied  to  all  orders.  Secondly,  an  account 
is  given  of  the  effect  of  a  velocity  slip  boundary  condition  on  the 
solution  of  the  Oseen  boundaiy  layer  equations  for  the  same  problem.  It 
is  found  that  an  additional  non -uniformity  in  the  solution  is  caused  by 
the  introduction  of  slip. 
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Drag  coefficient,  eq.  (44) 

7>ocal  skin  friction  coefficient,  eq.  (42) 

Constants  defined  in  eq.  (18) 

Constants  defined  in  eqs.  (23)  and  (30) 
functions  defined  in  eq.  (14) 

Functions  defined  xn  eqs.  (24),  (25),  (28)  and  (30) 
Functions  defined  in  eq.  (12) 

Functions  defined  in  eq.  (58) 

Functions  defined  in  eq.  (26) 

Constant  defined  in  eq.  (32)  ( *  ^as) 

Knuds en  number 

Constants  of  integration  (Appendix  II ) 

Symbolizes  terms  of  order  -  A/  ,  i.e.,  x *  0(x 

—  constant 
Pressure 


,  Reynolds  number  based 
UjaJL  ,  Reynolds  number  based 
7^  »  molecular  speed 

X  and  velocity  components 
Functions  defined  in  eq.  (46) 
Free  stream  velocity 
Cartesian  coordinates 
Constants  defined  in  eq.  (29) 
Constant  defined  in  eq.  (21) 


an  L 
on  x 
ratio 
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NOMENCLATURE  (continued) 
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Slip  parameter,  see  eq.  (36) 
■j-y-t  ,  auxiliary  variable 
tan"'  #/x 
Viscosity 


f*  /Q  t  Kinematic  viscosity 
Free  stream  kinematic  viscosity 
,  auxiliary  variable 


Density 

Free  stream  density 
Parabolic  coordinates,  eq.  (6) 
Stream  function,  eq.  (2) 

^  l**1?)  Functions  defined  in  eq.  (11) 
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Laplacian  operator 
Asymptotic  equality 


1 


1.0  INTRODUCTION 


Blasius  solved  the  incompressible  boundary  layer  problem  for  laminar 
flow  past  a  semi -infinite  flat  plate  at  zero  angle  of  attack  (references 
1  and  2)  and  his  result  describes  the  flow  quite  well  for  regions  where 
the  local  Reynolds  number,,  u^©  n  / V  ,  is  large.  Von  Karman  has 

proposed  (reference  3)  that  the  Blasius  solution  be  considered  the  first 
term  of  an  asymptotic  representation  of  a  solution  to  the  full  Navier- 
Stokes  equations;  calculation  of  the  second  term  should  extend  boundary 
layer  theory  to  lower  Reynolds  numbers.  Due  to  the  non-uniformities 
which  are  inherent  in  boundary  layer  theory,  any  such  representation 
will  contain  terms  which  express  these  non -uniformities;  in  particular, 
one  would  not  expect  a  simple  power  series  in  the  viscosity.  The  specific 
form  for  the  stream  function  expansion  suggested  by  von  Karman  has  a  line 
singularity  on  the  y-axis  (the  plate  is  assumed  to  occupy  the  positive 
x-axis  with  the  leading  edge  at  the  origin).  Alden  (reference  4)  improved 
this  situation  somewhat  by  modifying  von  Karman  *s  representation  in  such 
a  way  as  to  rotete  the  line  singularity  to  the  negative  x-axis.  Adopting 
von  Herman's  philosophy,  he  obtained  a  "first  correction"  to  the  Blasius 
solution.  At  first  glance,  Alden  's  iteration  scheme  seems  to  lead  to  an 
asymptotic  series  which  is  formally  a  solution  of  the  Navier-Stokes  problem. 
In  the  first  part  of  this  paper  it  will  be  shown  that  Alden 's  method  can¬ 
not  lead  to  such  a  solution  because  of  an  unavoidable  violation  of  the 
free  stream  boundary  conditions. 

In  the  case  of  a  gas  at  low  Reynolds  numbers  the  boundary  conditions 
at  a  solid  surface  must  be  modified  to  account  for  the  slipping  of  the  gas 
at  the  wall  (references  5  through  9).  To  support  the  conjecture  that  slip 
introduces  an  additional  non-uniformity  into  the  solution,  the  linearized 
boundary  layer  problem  (for  the  flat  plate)  with  a  slip  condition  is 
examined  in  detail  in  the  second  part  of  this  paper.  An  explicit  solution 
is  known  "and  this  is  compared  with  the  divergent  series  solution  obtained 
by  a  perturbation  in  the  "slip  parameter". 
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2.0  ANALYSIS 

2.1  Alden's  Method 

The  Navier-Stokes  equations  of  motion  and  continuity,  specialized 
for  the  steady  state  two  dimensional  flow  of  an  incompressible  fluid  are 
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where  x  and  ^  are  cartesian  coordinates,  u  and  x/-  are 
velocity  components,  is  the  Laplacian  operator,  and 

V  are,  respectively,  the  density,  the  pressure  and  the 
viscosity. 

A  stream  function,  ,  may  be  introduced  so  that 


(1) 


the  cartesian 
£  ,  p  and 

kinematic 


U  -  AJfc  (2) 

and  the  continuity  equation  is  automatically  satisfied.  If  the  momentun 
equations  are  cross  differentiated  and  subtracted,  to  eliminate  the  pressure 
teims,  one  obtains,  in  terms  of  the  stream  function 


If  =  <3) 

Physical  arguments  lead  one  to  suspect  that  the  boundary  layer 
equation  (reference  10) 

(4) 

may  describe  the  flow  quite  well  in  regions  wher  the  local  Reynolds 
nvmber,  ,  is  large,  (  u*>  being  the  free  stream  velocity) • 
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To  complete  this  stream  function  formulation  of  the  incompressible  boundary- 
layer  problem  for  laminar  flow  past  a  semi-infinite  flat  plate  at  zero 
angle  of  attack,  it  will  be  assumed  that  the  plate  occupies  the  positive 
x-axis  with  the  leading  edge  at  the  origin  and  it  will  be  required  that 

^  a  a  ^  ^  at  %  ~ 0  *  * >  0  ^ 

<X$  y  — >  o  <  9  <  T 

where  Y*  and  9  are  polar  coordinates  with  origin  at  the  leading  edge  of 
the  plate.  This  is  essentially  the  problem  solved  by  Blasius. 

In  terms  of  parabolic  coordinates, 

*  (<r  +  tx)a  *  (6) 

equation  (3)  becomes 
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Applying  the  analogous  boundary  layer  arguments  to  equation  (7)  one 
obtains  the  "parabolic  boundary  layer  equation", 
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Note  that  the  parabolic  boundary  layer  equation  is  not  the  same  as  the 
cartesian  boundary  layer  equation  expressed  in  parabolic  coordinates 
(see  reference  ll). 

H.  L.  Alden  (reference  4) >  following  a  suggestion  by  von  Karmen 
(reference  3),  attempted  to  improve  the  boundary  layer  solution  by 
developing  the  iteration  scheme  which  solves  equation  (l)  with  the 
Blasius  solution  as  the  first  term  in  an  expansion.  Using  his  scheme, 
Alden  found  two  terms  in  a  proposed  expansion.  The  first  term  is  the 
solution  of  the  parabolic  boundary  layer  problem  and  the  second  term 
may  be  thought  of  as  a  correction  term.  Since  there  is  no  obvious 
mathematical  reason  for  concluding  that  his  scheme  will  not  lead  to  an 
asymptotic  expansion  which  is  formally  a  solution  of  the  problem,  it 
is  worth  while  to  investigate  further.  Alden  introduces  the  auxiliary 
variables 

] 
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(9) 


and  rewrites  equation  (7)  as 


Alden  proposed  the  iteration  scheme  whereby  is  formally 

expressed  as 

f  (*>*,,*)  =  ?)  *  '>*¥'.<*>'})  *  v* %(«>'>)  *  -■  (11) 

and  this  is  to  be  substituted  into  equation  (10)  and  the  coefficients 
of  explicit  powers  of  V  set  equal  to  zero.  The  function,  (&>*))  > 

will  be  found  to  be  a  solution  of  the  parabolic  boundary  layer  problem 
and  this  scheme  may  be  thought  of  as  a  perturbation,  in  >>  ,  around  the 
parabolic  boundary  layer  solution.  If,  in  equation  (7),  V  is  set 
equal  to  zero,  the  order  of  the  equation  is  depressed.  One  therefore 
calls  this  perturbation  scheme  singular  and  expects  some  sort  of  singu¬ 
larity  in  the  solution  (references  12  and  13).  Furthermore,  the  ex¬ 
pansion  (11)  is  not  a  true  power  series  in  >>  since  the  coefficients 
%  l  r>  >?;  depend  implicitly  on  V  through  y  . 

Aldan's  scheme  produces  a  system  of  equations  for  the  ( <Tj  *7) 
and  for  each  equation  the  variables  may  be  separated  if  one  sets 
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(12) 


This  yields  a  system  of  ordinary  differential  equations  for  the  iff)  • 
Unfortunately  expression  (12)  implies  that,  for  n  £  2  , 

— ►  00  as  0“->  0  .  It  might  be  possible  however,  to  show 

that  the  expression  (11)  is  an  asymptotic  series,  uniform. in:  <r  and 
y  ,  in  some  region  excluding  an  arbitrarily  small  region  bounded  ty 
a  parabola  enclosing  the  stagnation  line  (  or  #  O  )  and  the  leading  edge. 
This  introduces  the  further  complication  of  whether  the  boundary  conditions 
which  remain ,  are  sufficient  for  the  problem.  The  region  mentioned  above 
excludes  the  origin,  which  is  the  most  likely  point  where  a  non -uniformity 
will  occur  in  the  solution. 

The  iteration  scheme  yields  the  following  equations  for  the  fj*) 

(i3) 

and  for  n  =  2,  3,  4,  .  .  . 
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where  is  dependent  on 

f  ,  -f,  ,  •  •  •  and  f„_,  ^ 

and  the  primes  indicate  differentiation  with  respect  to  J  . 

(  FiCf)  ,  F, t  ( t  and  ^4  [  ?)  are  exhibited  in  Appendix  I). 

The  boundary  conditions  (5)  require  the  following  behavior  for  the 
m)  ,  introduced  in  equations  (11)  and  (12), 
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From  equation  (13) 

-ft  +  -fj  f(W  '  ®  Cl(  a  constant.  . 


Alden  states  that  because  of  conditions  (16) 
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This  can  be  proved  by  making  the  highly  reasonable  assumption  that  the 
velocity  gradients  are  continuous.  Because  of  this  hypothesis,  ■£, 
and  its  first  two  derivatives  are  continuous  on  O  4  ^  .•  Now, 

rewrite  equation  (17)  as  follows. 
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Using  conditions  (16),  it  is  easy  to  show  that 
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Thus,  the  equation  for  ^  becomes 

C  +  f,f,’  *  0  (2°) 

This  equation  with  the  conditions  in  equations  (16)  possesses  a 
solution  (reference  14),  often  referred  to  as  the  Blasius  function, 
which  is,  together  with  all  of  its  derivatives,  continuous.  Furthennore, 
this  is  the  only  solution  with  continuous  derivatives  (reference  14)* 

The  Blasius  function  exhibits  the  following  asymptotic  behavior  as 
^  oc  (reference  15). 
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Using  equations  (6),  (9)  and  (12),  one  may  formally  rewrite  expression  (ll); 
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If  conditions  (16)  are  satisfied  and  each  fJV  is  continuous 
O  <  ^  <  *0  ,  then  Uoo  -fh  (%)  jz  ^  is  bounded  on  O  ^  f  <  •©  and 
expression  (22)  is  an  asymptotic  expansion. 

Unfortunately,  Alden  's  scheme  does  not  lead  to  a  solution  of  the 
Navier-Stokes  problem  since  )  cannot  be  made  to  satisfy  conditions 


This  will  be  shown  in  the  following  pages. 
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One  may  rewrite  equation  (14)  as 
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Integrating  again, 
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for  2,  3,  4,  . . . 
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Multiply  equation  (25)  by  /iH  f'jj)  end  obtain 


One  can  show  that 
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«*  M?)and  S-lt)  are  asymptdtically  finite,  i.e.,  for  each 
M  there  is  an  integer  ^  such  that  — *  0  as  J  . 
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Carrying  out  the  obvious  algebra,  one  can  show  (see  Appendix  II) 

that 


Expression  (32)  is  identical  with  Alden's  result  and  is  introduced  here 
to  show  that  K  0  .  In  fact  Alden  computed  K  =1.660  by  numerical 
integration.  Thus,  )  cannot  be  made  to  satisfy  condition  (16) 

and  consequently  the  proposed  expressions  (11)  and  (12)  cannot  be  made, 
even  formally,  to  satisfy  equations  (10)  and  (5).  It  may  be  remarked 
that  the  Alden  approach  with  a  slip  boundaiy  condition  is  subject  to  the 
same  difficulty. 


2.2  The  Flat  Plate  Boundaiy  Layer  Problem  With  Slip 

In  an  attempt  to  incorporate  rarefaction  effects  into  boundary  layer 
theory,  one  is  led  (reference  9) to  consider  a  slip  boundaiy  condition, 

,  where  ^  will  be  called  the  slip  parameter.  "J  is 

*3- 
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proportional  to  the  mean  free  path  length. 

The  flat  plate  boundary  layer  problem  with  slip  (i.e»,  with  a  slip 
boundary  condition)  has  not  yet  been  solved.  With  the  idea  of  gaining 
some  knowledge  of  the  slip  effect  in  general,  the  corresponding  linearized 
problem  for  incompressible  flow,  which  has  been  solved  explicitly  (refer¬ 
ence  7),  will  now  be  examined.  First,  the  problem  will  be  formulated  and 
the  solution  displayed.  Secondly,  a  series  solution  of  restricted  validity 
will  be  obtained  by  a  perturbation  scheme  in  ^  .  Lastly,  these  solutions 
will  be  compared  and  analyzed. 

The  linearized,  incompressible,  steady -state  boundary  layer  problem 
with  slip,  for  the  semi -infinite  flat  plate  at  zero  angle  of  attack  may 
be  stated 
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u  — >  u<*>  o 5  ,  x  >  o  (35) 

u  ”  T  y,  *  O  >  x  >  O  (36) 


This  is  reminiscent  of  the  Rayleigh  impulsive  plate  problem,  and  if 
condition  (35)  is  replaced  by 


U~Moo  at  x  ;  O  }  y.>0  |  (37) 

U.  bounded  in  the  entire  plane  i 

one  may  solve  this  new  problem  with  the  aid  of  Laplace  transforms  (refer¬ 
ence  7).  Fortunately,  the  solution  to  this  second  problem  satisfies  con¬ 
ditions  (35)  so  the  solution  to  the  original  problem  is  at  hand.  This 
solution  is: 
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If  one  defines  a  local  skin  friction  coefficient, 
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one  obtains,  for  this  problem, 


(40) 


(a) 


(42) 
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The  drag  coefficient  for  that  portion  of  the  plate  extending  from  the 
origin  to  a  distance  L  downstream  will  be  defined  (including  both  upper 
and  lower  plate  surfaces)  as 


(44) 


and  for  this  problem 


Equations  (38),  (43),  (45)  constitute  the  solution  of  the  originally 
posed  problem. 

An  attempt  to  solve  this  same  problem  by  a  formal  iteration  in  ~S  , 
commences  with  the  formal  series 
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This  yields  the  iteration  scheme 
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The  technique  employed  in  solving  the  problem  delineated  by  equation  (34) » 


under  the  conditions  (35)  and  (36)  may  be  employed  here  to  solve  equations 
(47).  The  no-slip  solution,  u6 3)is 


It  is  easy  to  see  that 


solves  equations  (48).  Thus,  formally, 


for  W  »  1,  2,  3,  .  .  .(50) 
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From  the  definition  (42)  and  equation  (51)  one  obtains,  formally 
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f j-rrf^r 


Integration  by  parts  and  the  use  of  induction  on  expression  (43)  yields 


2-  [3T 


ZHT(2p)('£f'j 


mo 


(53) 


Thus,  one  sees  that  the  formal  iteration  scheme  led  to  an  asymptotic 
(around  X  =  O  )  expansion  for  Cf  (which  is  uniform  in  x  for  x  bounded 
away  frcm  zero).  Expression  (52)  is  useless  at  x  =  O  and  cannot  help 
in  the  determination  o tCtf  it  can,  however,  in  a  series  of  unjustified 
steps,  be  summed  exponentially,  to  yield  expression  (43)  (see  Appendix  III). 

From  expression  (43 )p  one  learns  that 


c< 

_  7 

for  fixed  ~S  >  O 

x  — V  O 

l 

(54) 

S-*o 

ii 

?i?l 

for  fixed  X  *>0  J 

and  the  non-unifomity  of  in  and  x  at  T  =  X  =  ^  is  apparent. 
From  expression  (53), 


_z  r^f,  _  _l 


(55) 


Integration  ty  parts  and  the  use  of  induction  on  expression  (45)  yields 


From  equations  (55)  and  (56)  it  is  evident  that  a  slip  correction 

1 

to  the  local  skin  friction  of  order  i  does  not  imply  that  the  slip 
correction  to  the  drag  is  also  of  order  >  since  the  non-uniformity 
is  sufficiently  strong  so  that  the  actual  correction  is  of  order  3  • 


3.0  CONCLUSIONS 


1.  The  scheme,  which  Alden  used  to  correct  the  Blasius  solution,  cannot 
lead  to  a  solution,  even  formally,  of  the  incompressible  Navier- 
Stokes  semi-infinite  flat  piste  problem. 

2.  For  the  incompressible  linearized  flat  plate  problem  the  inclusion 
of  a  slip  boundary  condition  introduces  a  correction  to  the  local 
skin  friction  of  order  'S*  but,  due  to  a  non-uniformity  in  the 
solution,  corrects  the  drag  coefficient  to  order  . 


APPENDIX  I 


r3(y)  =  4sr-rr-  ru“ - iuc + ra: 


MV)  = 


-  xc  -  w  -  <*  -  nr  -  -  m 

->H‘i  -  tfX  +  nC -  «£* 

+*u"  -  r-c 


w)  -  -m;  u.r  -xc  -fX  +mx~hX 

~6°x  ->nr-nx  -m;  -m** 
-xa  -nr +aa'+rum-*a‘-iu'{ 

*n:r^nx  -hu'  +*xi  -*m 

-  *  T-Cf'  -  ^  if  -  W/  +  *  X  -  fC  -  st{ 
+  *X'-XA‘ 


17 


APPENDIX  II 

From  Appendix  I  and  equation  (21)  F^  0  «  Introducing  the  integration 
constants,  ,  one  obtains 

> — ■  o 

f;,  —  + 


r  v— ^ 

(ca,  +  k3,)C  j-^/s)  * 

C„  +  Kaa 

a.  - 

i 

(ii-D 

(ca,  ♦  Kai)  d^i)  + 

(Caa.  +  K„)(5 

~/3/a)  -t  K*3 

^3 

f»  - 

»  Cji^KaiT  ^  “/^/a) 

a  L  2. 

+  o(,n]  - 

f-  C»x  +  Kaa 

2. 

+  !^[cs-fo)"  ■ 

+  ofr') 

Because 

of  conditions  (16) }  one  i 

must  put 

* 

'-a » 

-  K„ 

(11 -2) 

Using  equations  (14) »  (24)  > 

(25)  and  (28), 

one  obtains 

f*  1 

^  K;1  +  Kal[ 

-/Va)'1  .  ("5  - 
"2  < 

]  +  o(v') 

4  -- 

3(5  -(3/a)'*j 

♦  °(V‘) 

c  - 

*>,  [(1  -/Va)'5  +  ^  A*)"*  ]  + 

o(r) 

(II -3) 

h  ' 

1  * 

-  °Tr') 

-m 

»a 

-  K„[l 3(f-/Vz)'f  + 

to  (5  -A/af 

+  o(r’) 
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Appendix  I  and  equations  (21)  and  (II-3)  yield 

f\  =  -io(C3a  +  K32)  +  6fi  Ka3  (5  -/Va)'a  -+  cCS'5) 

Proceeding  in  the  same  way  as  with  Fx  , 

F3,  =  -'OCCa  +  C„  -  Kj, 

->  0(5“) 

F,i  =  *«)($- /3/a)'  +  (C„  +  Ki,)(5-fih) 

-6fiKn  4*  (?-/!/»)  +0(5") 

Fn  -  -*  (C*  +  Ktt)(j-fl/a)r[l  t  £*r*)] 

+  4r„  +  k„  0  ( S'1)] 

-8/SK.»(1-M)’A^S-/3/4)f  I  *  0(Y*)]  + 

fMC  =  ii<5'W  +  °fr‘) 


(II -4) 


(II-5) 


°(r3) 
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^  +  CS±iS>[h-  0(5-)] 

'  i  A  K„  ( Y-fl/tf'  A.  (r-M)fi  +  0(5- jj 

•*  o(r) 

£  =  -  i(43  +  *0  ( 5  -fl/*)'[  I +0(  5-)] 

+  (Y-fi/i)  [ i  *  0(5-)] 

-  (s  a  ( s-flA)  [  i +  <>(rv] 

+  0(,) 

Because  of  conditions  (16), 

^ai  r  “  Kaj 
^Ji  “  “^*1 

£  *-/JKnAft-Afc)  -f  a,(5-A0°(V*) 

£  =  -/JK„(5-/o/l)"  -  ofV1)  +A  (5-/»A)o(s-») 

{/'  =  /J  K„  (5  -/Va)'a  -  0(5°)  *  A  )o(5‘») 
f“  - -2/3K„/S'/3/»)'5  *  «(V4) 

{“  =  6flKt,  (5 -M)'4  *  0(F*r)  ♦  A,  (  5  -A'a)o(s-) 


(II-6) 


(II -7) 


(II -8) 
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fr  =  +o0)  +  ^  ’) 

fn  =  °(0 

+  0(S*)  +  Ai  (5-fi/t)  °(S) 

Ut  *  ilte-fl/i)*  -‘fStS-ft/iY  *  n 
F„  -  m  fiKntS-  /Vo)’  MS-  ^O' 

■  +[j^fs-fa)Yo(r)  +  oiv) 

*  A  fJ-M)  °(V) 

fv4  *  i 

+  (A'(j-m)TW)  +°(r‘) 

■*  Jt*(S-fl/t)o(5-*) 
f»  =  fyS  K„  f  ?- A4)’-<2~  /  ?-/»/>) 

A. /*-««]**'>  +  °/t*) 

+  A.(V0/SWS) 

[h(s)T  -  ikH-Arf  +  ^r") 


(II -9) 
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APPENDIX  III 

Proceeding  from  equation  (52)  in  a  purely  formal  manner: 


The  expression 


satisfies 

~  2*1  +  2  *  O 

1(6)  -  JrT 

Therefore 


(III-1) 


(III -2) 


(III-3) 


(HI-4) 


X  (k)  -  */rF  (  X')  (  M 


(III -5) 
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and  finally, 

4  = 

This  method  of  “exponentially  summing*  a  divergent  series  is  discussed 
ty  Borel  in  reference  16. 


m)~* 


(III-6) 
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